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Abstract. We present an investigation of the effects of asphericity on the estimates of total mass and gas mass 
fraction in galaxy clusters from X-ray observations. We model the aspherical shape of galaxy clusters by a triaxial 
model and compare the true total mass and the true total gas mass fraction with the corresponding quantities 
obtained with the assumption of spherical symmetry. In the triaxial model we allow the extent along the line 
of sight to vary in order to describe elongated and compressed cluster shapes. Using a sample of 10 ROSAT 
clusters and a recent CHANDRA observation we find the following results. For prolate or oblate shapes the 
difference between triaxial and spherical model both in the mass and in the gas mass fraction are negligible 
(less than 3%). For more aspherical shapes the total mass is underestimated (overestimated) in the centre, if the 
cluster is compressed (elongated). The gas mass fraction is overestimated for compressed clusters and slightly 
underestimated for elongated clusters. Comparing X-ray masses with gravitational lensing estimates, we find that 
elongations along the line of sight can resolve discrepancies of masses determined by the two different methods of 
up to ~ 30%. The combination of Sunyaev-Zel’dovich and X-ray observations is useful to measure the elongation 
of the cluster along the line of sight. As an application, we estimate the elongation of the cluster CL0016-I-16 with 
two different approaches, Sunyaev-Zel’dovich measurements and comparison of weak lensing and X-ray masses, 
and find reasonable agreement. 
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1. Introduction 

Clusters of galaxies are the largest gravitationally bound 
aggregates of matter in the universe and, for many 
applications, they can be regarded as being representative 
for the universe as a whole. In particular the ratio of 
baryonic to total matter in clusters can be assumed to 
be representative of the universe, because of the large 
volume and the fact that clusters are closed systems. 
Combined with big bang nucleosynthesis calculations 
and observed light-element abundances, this ratio can be 
used to constrain the cosmological density parameter Qm- 
Optical velocity dispersion measurements and gravita¬ 
tional lensing are independent methods used to constrain 
the total mass in galaxy clusters while X-ray observations 
and recently the Sunyaev-Zel’dovich (SZ) effect (Grego et 
al. 2001) are used to estimate not only total masses but 
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also gas mass fractions in clusters. The X-ray method is 
based on the assumptions that the X-ray emitting gas is 
an ideal gas in hydrostatic equilibrium and the total mass 
and gas mass fraction are estimated through the X-ray 
surface brightness distribution and the gas temperature. 
In this context, the /3-model (Cavaliere & Fusco-Femiano 
1976) is widely used in X-ray astronomy to parametrise 
the gas density profile in clusters of galaxies by fitting 
their surface brightness. It is usually used under the 
assumption of spherical symmetry. If the images show 
that the cluster emission is smooth and the isophotes 
show an elliptical shape, the modelling is improved when 
an elliptical /3-model is used: Fabricant et al. (1984) 
studied the non-spherical shape of A2256, McMillan et 
al. (1989) the morphology of 49 Abell clusters and the 
elliptical shape of C100I6-I-I6 was analysed by Neumann 
& Bohringer (1997) and Hughes & Birkinshaw (1998). 
Buote & Canizares (1996) studied the elliptical shapes of 
5 Abell clusters concluding that, for oblate and prolate 
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shapes, the estimates of the gas masses are insensitive to 
the ellipticities of the X-ray isophotes. 

The aim of this paper is to investigate to which degree 
the assumption of spherical symmetry in the estimate 
of total mass and gas mass fraction in galaxy clusters 
showing elliptical X-ray emission is accurate. We select 
10 galaxy clusters showing a smooth, elliptical emission 
in the ROSAT archive and we model the intracluster 
medium (ICM) with a triaxial /3-model and compare 
the estimates with those of a spherical /3-model. In 
particular, a triaxial model allows the description of 
observations of clusters which may be compressed or 
elongated along the line of sight, with the latter being 
probably more frequent, since elongated clusters are se¬ 
lected preferentially due to their higher contrast over the 
background emission. Since our investigation is focused 
on a geometrical aspect common to many clusters, we 
do not present an exhaustive discussion on the results 
for each single cluster. Instead, we point out the general 
features common to each sample cluster. In Sect. we 
describe our model. In Sect. ||we present the sample and 
the determination of the morphology and the relevant 
best fit parameters necessary for the estimate of total 
masses and gas mass fractions, which are discussed in 
Sect. ^ and ||, respectively. We discuss the influence of 
asphericity on the comparison of X-ray and weak lensing 
mass estimates in Sect. || and investigate further possible 
constrains on the cluster shape by means of the SZ effect 
in Sect. |^. We present a summary and conclusions in 
Sect. 1^ 

If not stated explicitly, Hq = 50kms“^ Mpc“^ and 
go = 0.5 are assumed. Icr uncertainties for the estimated 
parameters are used throughout this paper. 


2. The model 


We extend the classical, spherical /3-model to a more gen¬ 
eral triaxial, ellipsoidal /3-model: 


/ x‘^ 


- 3 / 3/2 


( 1 ) 


for the electron density in the ICM, where a, b and c are 
the core radii in the x, y and z directions respectively. 
Assuming isothermality of the ICM, the X-ray surface 
brightness, Sy^, is proportional to the integral of the square 
of the electron density along the line of sight. Taking the 
latter to be along the z axis we obtain: 


S'x(x, y) (xcnl^x + ^ + 


- 30 + 1/2 


( 2 ) 


Thus, the core radii a and b can be determined directly 
from X-ray images, while some assumptions on the core 
radius along the line of sight c must be made. On the other 
hand, if also a possible inclination of the principal axes of 
the isodensity surfaces described by this triaxial model is 


Table 1. The values for the basic physical parameters 
compiled from the literature: the redshift z and the gas 
temperature Tgas. References: [I] Hughes & Birkinshaw 
(1998), [2] Matsumoto et al. (2000), [3] Allen & Fabian 
(1998), [4] Markevitch et al. (1998), [5] Ettori & Fabian 
(1999). 


Cluster 

z 

2gas 



(keV) 

C10016-tl6 

0.5455 

7.557;j'i^8[l] 

A 478 

0.0881 

6.40t°;“[2] 

A 1795 

0.0631 

5.68t°:lj[2] 

A 1068 

0.1386 

5.5l(,;^[3] 

A 1413 

0.1427 

8-5la8[3] 

A 2390 

0.231 


A 2199 

0.0298 

4.221^6 [2] 

A 2029 

0.0765 

8.47i°:«[3] 

A 2597 

0.0852 

3.6«J[4] 

Hydra A 

0.0522 

3.7lt“:ll[2] 


included, only a combination of the three core radii can 
be determined from the X-ray images. The investigation of 
this possibility, besides being involved, does not substan¬ 
tially improve the modelling and we do not expect that 
it will change our conclusions. We thus model the X-ray 
emission with a spherical model and a triaxial model with 
no inclination angle to the line of sight. For the spherical 
model a — b = c = in Eq. (|^ and for the triaxial 
model we choose the major core radius a = r®", the minor 
core radius b = e x r®" and the core radius along the line 
of sight c = i X r®^b As noticed by Grego et al. (2000) 
for A 370, the total mass density computed under the 
assumption of hydrostatic equilibrium and isothermality 
of the intracluster gas can be unphysical for very ellipti¬ 
cal gas density distributions. Consequently, only a specific 
range for the parameter i, which depends on the cluster 
properties, is allowed in the modelling of ellipsoidal ICM 
distributions with a triaxial, ellipsoidal /3-model. 

3. Sample and morphological analysis 

The clusters are selected from the ROSAT archive. Since 
our goals require a sample of clusters that are regular and 
show a smooth, elliptical surface brightness, we exclude 
all the clusters with bimodal or strongly irregular mor¬ 
phology and those with a high spherical symmetry. In 
order to have a robust estimate of the parameters that 
constrain the morphology, we select only those clusters 
for which both ROSAT High Resolution Imager (HRI) 
and Position Sensitive Proportional Counter (PSPC) ob¬ 
servations are available. Because of the lack of detailed 
temperature maps and for simplicity, we assume that the 
gas in each cluster is isothermal. The 10 selected clusters 
with the additional information needed to deproject the 
surface brightness are listed in Table |^. For both HRI and 
PSPC, we prepare exposure corrected images and remove 
point sources embedded in the cluster emission to prevent 
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Table 2. The estimated morphological parameters: the 
position angle (N over E) and the minor to major axis ratio 
e and the ranges of allowed elongations (imin < * < *max)- 


Cluster 

Position 

angle(deg) 

e 

^min 

^max 

C10016-tl6 

47.1 ± 7.7 

0.85 ± 0.06 

0.64 

1.80 

A 478 

47.4 ± 7.4 

0.84 ±0.06 

0.65 

1.56 

A 1795 

10.2 ± 6.4 

0.83 ±0.06 

0.64 

1.48 

A 1068 

-47.0 ± 4.9 

0.79 ±0.04 

0.62 

1.27 

A 1413 

-5.1 ±4.3 

0.72 ± 0.03 

0.59 

1.05 

A 2390 

-60.0 ±8.1 

0.79 ±0.07 

0.62 

1.31 

A 2199 

32.9 ± 9.2 

0.86 ± 0.04 

0.66 

1.68 

A 2029 

16.1 ± 7.4 

0.80 ±0.06 

0.63 

1.35 

A 2597 

-35.3 ±5.1 

0.79 ±0.06 

0.62 

1.30 

Hydra A 

-32.9 ± 7.1 

0.87 ±0.05 

0.66 

1.75 


contamination on measurements of the ellipticity of the 
surface brightness. The determination of the isophotes el¬ 
lipticity and position angle is performed with the MIDAS 
routine FIT/ELL3, an iterative least-squares method in 
which the isophotes are free to translate position, and 
change ellipticity and orientation (Bender & Mollenhof 
1987). From each image processed with this algorithm we 
obtain a set of ellipses with given position angles, half mi¬ 
nor and major axes and centre coordinates. Even though 
our model implies that the isophotes are similar and con¬ 
centric, such a set of parameters allows us to identify the 
isophotes which are suitable to describe the shape of the 
cluster emission. In the PSPC data we notice that for all 
the clusters the minor to major axis ratio e tends to unity 
as we approach the centre. We identify this effect as a 
distortion due to the point spread function (PSF) of the 
PSPC and thus discard the isophotes within the central 
~ 25 arcseconds region. The X-ray emission can be traced 
up to Rout, which is estimated from the surface bright¬ 
ness profile and listed in Table 5. For each cluster we also 
can identify a distance from the centre beyond which the 
isophote parameters jump to arbitrary values. As this dis¬ 
tance is approximately equal to Rout, we also discard the 
isophotes beyond it. We then compute the mean values of 
the minor to major axis ratio e, position angle and centre 
position from the parameters of the remaining isophotes. 
For the HRI images we also notice the same central fea¬ 
ture, but due to the smaller held of view, we are not able 
to estimate the ellipses parameters for the whole cluster 
emission. In any case we hnd that the mean values for the 
minor to major axis ratio e, the position angle and centre 
computed from HRI and PSPC images agree within the 
errors when evaluated over the same radial range. For the 
analysis we use the PSPC images because of the PSPC’s 
superior sensitivity and larger held of view. The relevant 
morphological parameters (the position angle and the mi¬ 
nor to major axis ratio e) from the image analysis are 
listed in Table |^. Since we are interested in the differences 
between spherical and ellipsoidal geometries in the esti¬ 
mate of total masses and gas mass fractions, we extract 


Table 3. The best ht parameters for the circular prohle: 


the central surface brightness the slope and the 

circular core radius r®'®®. 

Cluster 

c*circ 

/IQ— 5 PSPC cts \ 

^circ 

circ 
' c 

(arcsec), (kpc) 

C10016±16 

0.91 ±0.05 

0.81 ±0.03 

51.51 ±2.94, 379 

A 478 

12.35 ±0.36 

0.60 ±0.01 

39.91 ±0.98, 88 

A 1795 

12.94 ±0.24 

0.60 ±0.01 

51.76 ±0.83, 85 

A 1068 

12.27 ±0.73 

0.68 ±0.01 

25.73 ± 1.30, 83 

A 1413 

3.51 ±0.20 

0.66 ±0.01 

53.74 ±2.88, 176 

A 2390 

4.52 ±0.29 

0.65 ±0.01 

41.09 ±2.33, 192 

A 2199 

6.70 ±0.14 

0.57 ±0.01 

84.55 ±1.58, 70 

A 2029 

17.92 ±0.46 

0.57 ±0.01 

40.66 ±0.90, 79 

A 2597 

16.85 ±0.73 

0.65 ±0.01 

30.04 ± 1.09, 64 

Hydra A 

15.15 ±0.54 

0.57 ±0.01 

35.14 ±1.00, 49 


Table 4. The best ht parameters for the elliptical prohle: 
the central surface brightness the slope and the 

core radius along the major axis r®'^ 

Cluster 

-- 

*^0 

/IQ—5 PSPC cts \ 


- jiin - 

' C 

(arcsec), (kpc) 

C10016±16 

0.89 ± 0.04 

0.82 ± 0.03 

57.11 ± 3.24, 420 

A 478 

12.20 ±0.34 

0.59 ±0.01 

43.02 ±1.02, 95 

A 1795 

12.30 ±0.25 

0.60 ±0.01 

58.77 ±1.04, 97 

A 1068 

12.12 ±0.68 

0.68 ±0.01 

28.76 ±1.37, 92 

A 1413 

3.46 ±0.23 

0.64 ±0.01 

58.85 ±3.65, 193 

A 2390 

3.99 ±0.25 

0.65 ±0.01 

48.67 ± 2.84, 227 

A 2199 

6.73 ±0.14 

0.57 ±0.01 

89.31 ±1.66, 74 

A 2029 

18.44 ±0.53 

0.56 ±0.01 

42.84 ±1.06, 84 

A 2597 

16.80 ±0.81 

0.64 ±0.01 

32.72 ±1.31, 70 

Hydra A 

15.06 ±0.51 

0.57 ±0.01 

37.89 ±1.01, 53 


two surface brightness prohles from the PSPC images, re¬ 
ferred to as circular and elliptical prohles. The circular 
prohles are determined by computing the counts within 
circular annuli spaced by 5 arcsec, while for the elliptical 
prohles the bins are similar and concentric ellipses with 
the minor to major axis ratios e and position angles listed 
in Table 2 and spaced by 5 arcsec in the direction of the 
major axis. We parametrise the circular and elliptical pro¬ 
hles with: 

Sx{r)=Si(l+[-) ) +BG\ (3) 

where j = circ, ell, labels the circular and elliptical mod¬ 
els, respectively. For the elliptical prohle, the distance r 
from the centre is measured along the major axis and 
thus the elliptical core radius r®" lies along the major 
axis. The best ht parameters are listed in Tables ^ and 
We notice that the central surface brightness of the 
elliptical prohle is in general smaller than that of the 
circular prohle and that the values for j3 are almost the 
same for both models. Obviously the elliptical core radius 
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Table 5. rsoo, the outer radius i?out and the total masses 
and gas mass fractions for the spherical model for R = 

^ 500 - 


Cluster 

rsoo 

(Mpc) 

-Rout 

(rsoo) 

MlotirU 

(WMo) 


C10016-F16 

1.051 

2.1 

6.23 ±0.93 

0.268 ± 0.027 

A 478 

1.492 

1.5 

6.22 ± 0.34 

0.256 ±0.019 

A 1795 

1.456 

1.2 

5.39 ±0.20 

0.217 ±0.015 

A 1068 

1.384 

1.2 

5.69 ± 1.55 

0.131 ±0.018 

A 1413 

1.677 

1.9 

10.23 ± 1.73 

0.157 ±0.013 

A 2390 

1.731 

1.9 

14.05 ± 2.18 

0.186 ±0.015 

A 2199 

1.288 

1.3 

3.39 ±0.11 

0.179 ±0.013 

A 2029 

1.713 

1.1 

9.11 ± 0.61 

0.212 ± 0.019 

A 2597 

1.172 

1.4 

2.99 ± 0.21 

0.183 ±0.011 

Hydra A 

1.170 

1.3 

2.71 ± 0.15 

0.207 ±0.017 


is always larger than 


4. Results for the total mass 



0 0.5 1 1.5 2 

R [rsool 


Fig. 1. The relative errors for the total mass estimates 
for A2390 plotted for the 4 triaxial models: compressed 
along the line of sight (z = Zmin) (dotted line), prolate 
(z = e) (dashed-dotted line), oblate (z = 1) (dashed line) 
and elongated (z = imax) (solid line) shapes. Positive rela¬ 
tive errors imply overestimations of the total mass if spher¬ 
ical symmetry is assumed: this is the case if the cluster 
is elongated. Underestimations are found for compressed 
clusters. 


Assuming hydrostatic equilibrium and isothermality of the 
intracluster gas, the total mass density ptot of a cluster can 
be computed from the /3-model with the best fit parame¬ 
ters listed in Tables ^ and 


where G is the gravitational constant, /cb the Boltzmann 
constant, Pgas the gas mass density, Tgas its temperature 
and pzTZp is the mean particle mass of the gas (we as¬ 
sume p = 0.61). In order to compare the total masses for 
spherical and ellipsoidal ICM shapes, we compute them 
by integrating Eq. (^) within a sphere of radius R. We 
choose to integrate the total density within a sphere for 
both models because, besides being the natural choice for 
the spherical modelling, it allows the comparison of masses 
contained within the same volume. We quantify the dif¬ 
ference between the two models by defining the relative 
error: 


EU{R) = 


Mll^iR) - MUR) 


MUR) 


( 5 ) 


where mU{R) and Ml^^{R) are the total masses within 
the distance R from the centre for the spherical model 
and a triaxial model with a core radius along the line of 
sight c = z X r®", respectively. Since the ellipsoidal mod¬ 
elling should be more accurate than the spherical, positive 
(negative) values of are going to be characterised as 
over- (under-) estimations. For spherical symmetry total 
masses evaluated from the best fit parameters of the cir¬ 
cular analysis for R = rsoo (the radius which encompasses 
a volume that has a total mass density 500 times the criti¬ 
cal density for closure) are listed in Table ||. In the case of 
ellipsoidal shapes some constrains on the core radius along 


the line of sight, c = z x r®", must be taken into account. 
As outlined in Sect. the total mass density inferred by 
means of Eq. (^ can be unphysical for certain ellipsoidal 
shapes of the ICM. In fact, very elliptical gas density dis¬ 
tributions would imply regions with negative total mass 
density. As we model the cluster within a sphere with ra¬ 
dius i?out, we clearly require the total mass density to be 
be positive inside such a sphere. Consequently, we obtain 
a range for the parameter z (zmin < * < *max) for each 
sample cluster. These values are listed in Table Notice 
that very ellipsoidal shapes are not allowed. We emphasise 
that the computed range Zmin < * < *max implies negative 
values for the total density only at distances abundantly 
larger than i?out, where anyhow we do no longer expect 
the adopted triaxial model to be a reasonable approxima¬ 
tion. Indeed, at very large radii the shape of the cluster is 
dominated by the infall of galaxies and groups of galax¬ 
ies. Therefore, the cluster shape there has certainly no 
similarity to the shape of the inner region, and hence the 
problem of a negative total mass density does not arise. 
We present the results of our investigation of the relative 
errors E\j^ for: the two axisymmetric ellipsoidal geome¬ 
tries; oblate (c = r®", z = 1 in Eq. (^)) and prolate (z = e) 
spheroids, a shape for which the distribution is the most 
compressed along the line of sight (z = Zmin) and a shape 
that describes the most elongated shape (z = Zmax)- For 
simplicity, we show in Fig. the relative errors i3J^(i?) 
for these 4 ellipsoidal models only for the cluster A2390, 
since the trend of these functions is qualitatively the same 
for all the sample clusters. We notice that the errors for 
the total masses in Table 5 include the uncertainties on 
the gas temperatures, which are quite large. Instead, if 
this is not included, the errors are of the order of 2 — 3%. 
Consequently, errors of this order which are solely due to 
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Table 6. The relative error for the total masses for R = 

’" 500 - 


Cluster 

ZT’^min 

-^M 

(%) 

Rm 

(%) 

Rm 

(%) 

ZT’^max 

(%) 

C10016-tl6 

-3.44 

-1.36 

-0.08 

4.72 

A 478 

0.84 

0.89 

0.93 

1.04 

A 1795 

-0.65 

-0.60 

-0.55 

-0.45 

A 1068 

0.73 

0.78 

0.83 

0.89 

A 1413 

3.73 

3.84 

4.06 

4.09 

A 2390 

0.32 

0.52 

0.73 

1.00 

A 2199 

0.85 

0.90 

0.93 

1.04 

A 2029 

2.05 

2.08 

2.10 

2.15 

A 2597 

1.51 

1.55 

1.59 

1.65 

Hydra A 

-0.20 

-0.17 

-0.16 

-0.08 


Table 7. The relative error for the total masses for R = 

„chc 

I r' 


Cluster 

{%) 

Rm 

(%) 

Rm 

(%) 

IT'^max 

-^M 

(%) 

C10016-tl6 

-11.62 

-1.98 

3.52 

20.23 

A 478 

-11.64 

-3.34 

2.29 

15.05 

A 1795 

-9.87 

-0.93 

5.65 

17.79 

A 1068 

-12.19 

-4.65 

3.16 

10.28 

A 1413 

-14.24 

-8.45 

0.96 

2.31 

A 2390 

-7.68 

1.26 

9.73 

18.69 

A 2199 

-12.18 

-3.83 

1.03 

15.40 

A 2029 

-14.25 

-7.07 

-0.55 

7.70 

A 2597 

-13.20 

-5.70 

1.61 

9.15 

Hydra A 

-8.13 

-1.67 

1.77 

14.08 


the non-spherical modelling of the cluster are irrelevant. 
We thus emphasise that only the relative errors 
larger than ~ 3% are significant. 

Since most of the mass is contained in the core, we find 
that at large radii, at r^Qo for instance, the relative errors 
if are almost independent on i and have values of a few 
per cent. This means that at large distances from the cen¬ 
tre the possible compression or elongation along the line of 
sight is not important and that the assumption of spherical 
symmetry yields very reliable results. The relative errors 
for R = rsoo for the 4 models are listed in Table As R 
decreases the relative errors in the mass estimations get 
larger (see Fig. for shapes elongated along the line of 
sight we always find overestimations which increase with 
the value of the core radius along the line of sight. The 
reason is that for strong elongations a large portion of the 
core is excluded from the sphere within which the mass is 
estimated. A similar explanation holds for the underesti¬ 
mation found for compressed shapes. In Fig. |^, we show 
the relative errors at R = for the 10 sample clusters as 
a function of z. At i? = we find: < >s= —12 %, 

< >s= -4%, < RIj >s= 3%, < >s= 13%, 

where < >s means that we average over the 10 sample 
clusters. As the masses are usually estimated at i? = rsoo 
or larger, these errors usually do not affect the mass de- 
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Fig. 2. The relative errors for the total mass R^j as a 
function of i for the 4 models (zmin, e, 1 and Zmax, which 
are separated by the vertical lines) at for the 10 sam¬ 
ple clusters. For clusters elongated along the line of sight 
we find overestimations, while for compressed clusters the 
total mass is underestimated. 

termination. 

We also perform the same analysis on a recent CHANDRA 
observation of the galaxy cluster RBS797 (Schindler et al. 
2001), which shows an elliptical emission with a minor 
to major axis ratio of 0.77. This cluster is an excellent 
object for this analysis because the minor to major axis 
ratio e and the position angle do hardly change with ra¬ 
dius. For this cluster we find MjQ^^(r 5 oo) = 6.5 x IO^'^Mq 
and /|P^(r 5 oo) = 0.17 and we arrive at the same conclu¬ 
sions on the relative errors for the total mass and gas mass 
fraction obtained from our ROSAT sample. For instance, 
comparing these values with the estimates from the 4 tri- 
axial models discussed in this section, we find that the 
relative error for the total mass within rsoo is always less 
than 1%. 

5. Results for the gas mass fraction 

We investigate the same 4 triaxial models as in Sect. ^ 
with respect to the gas mass fraction. For the computation 
of the gas mass fractions for the 10 sample clusters we 
deproject the X-ray emission with the /3-model in order 
to determine the gas density pgas and then integrate it 
within a sphere of radius R. We thus obtain the gas mass 
fractions fg^{R) for the spherical geometry, which are 
listed in Table || for R = rsoo, and the gas mass fractions 
/gas(^ for the triaxial models, where i is defined as in 
Sect. y. Similarly, we define the relative error: 

/fP^R) - ftasiR) 

■ — Turn —■ ‘' 

For the reason discussed in Sect. |^, only relative errors 
for the gas mass fractions larger than ~ 3% are of 
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R [Tsoo] 

Fig. 3. The relative errors for the gas mass fraction es¬ 
timates for A2390 plotted for the 4 triaxial models: com¬ 
pressed along the line of sight (z = Zmin) (dotted line), pro¬ 
late (z = e) (dashed-dotted line), oblate (z = 1) (dashed 
line) and elongated (z = Zmax) (solid line) shapes. For com¬ 
pressed clusters the gas mass fraction is overestimated. 


significance. The general trend of the relative errors 
Ef{R) is the same for all the sample clusters and can be 
seen in Fig. where the relative errors for the cluster 
A2390 are shown. 

Since the gas mass fraction fg^s{R) is the ratio of the gas 
mass Mgg^{R) to the total mass Ml^^{R), the discussion 
of the relative errors Ef{R) involves many quantities. 
We notice that, since the central surface brightness is 
proportional to c x PgasOi the central gas density is 
larger for the compressed shapes than for elongated. 
This factor is usually dominant for small R, leading to 
gas mass overestimations for elongated ellipsoids and 
underestimations for compressed shapes. In Table we 
list the gas mass fraction relative errors for the 4 models 
at i? = . We find small over- and underestimations for 

oblate and prolate shapes for small R: < \Ef\ >s= 2% 
and < \Ef\ >s= 3% at R = For elongated and 

compressed shapes we find larger but in general still 
negligible errors: at i? = we get < >s= 6% 

and < >s= 5 %. 

At large distances the core radius along the line of 
sight is dominant: small core radii imply steep profiles 
and consequently less gas mass. In Fig. we show the 
relative errors at i? = rsoo for the 10 sample clusters 
as a function of z. We see that at large distances from 
the cluster center, at rsoo or larger, where the total 
masses for a spherical and an elliptical model are almost 
the same (see Sect. the compressed shapes imply an 
overestimation of the gas mass fraction for the majority 
of the clusters. Although this overestimation is small, the 
more compressed the shape, the larger the overestimation 
(see Fig. ^). At i? = rsoo we find the average values: 
< >s= 3%, and < Fif*"** >s= -2.8% at i? = rgoo, 

showing small underestimations of the gas mass fractions 
for elongated shapes (see also Table ^). For oblate (z = 1) 
and prolate (z = e) shapes we find that at i? = rsoo 


Table 8. Relative 
R = 

errors 

for the 

gas mass fractions 

Cluster 

^^min 

(%) 

Ef 

(%) 

Ej 

(%) 

^^max 

(%) 

C10016-tl6 

7.54 

1.37 

-0.29 

4.94 

A 478 

4.91 

1.90 

1.38 

6.19 

A 1795 

3.96 

0.64 

-0.04 

3.70 

A 1068 

6.78 

2.94 

1.57 

2.94 

A 1413 

10.74 

7.73 

6.57 

6.81 

A 2390 

5.31 

1.43 

0.24 

2.05 

A 2199 

4.53 

1.83 

1.52 

8.06 

A 2029 

6.40 

3.88 

3.56 

6.48 

A 2597 

7.44 

3.93 

2.94 

4.82 

Hydra A 

2.30 

0.40 

0.40 

7.92 

0.6 

0.8 

1. 

1.2 

1.4 1.6 1.8 


10 

7.5 

5 


LU 0 
-2.5 


0.6 0.8 1. 1.2 1.4 1.6 1.8 

i 

Q016 + 16A478 A1795 A1068 A1413 A2390 A2199 A2029 A2597 HydraA 
OOOx □A*0 + ■ 


Fig. 4. The relative error for the gas mass fraction El as a 
function of z for the 4 models (zmin, e, 1 and Zmax) at rsoo 
for the 10 sample clusters. The 4 elongations are separated 
by 3 vertical lines. 

both underestimates and overestimates are present and 
that they are small (see Table ^. Averaging the absolute 
values of the relative errors we find: < \El\ >s= 1.7% 
and < |i7f I >s= 2.6%. We conclude that these errors 
usually do not substantially affect the gas mass fraction 
determination at large radii. 


6. Asphericity and comparisons with gravitational 
tensing 

Since the X-ray method and gravitational leasing are inde¬ 
pendent methods used to determine the dark matter dis¬ 
tribution in cluster of galaxies, their comparison is very 
useful. First of all, the morphology of the dark matter 
distribution can be compared: Smail et al. (1995) found 
good agreement for the position angle and ellipticity of 
C10016-I-16 determined with weak leasing and X-ray ob¬ 
servations. Furthermore, Miralda-Escude & Babul (1995), 
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Table 9. Relative errors for the gas mass fractions at 

R = ''’ 500 - 


Cluster 

^*min 

(%) 

Ef 

(%) 

E} 

(%) 

(%) 

C10016-F16 

8.57 

2.52 

0.30 

-0.13 

A 478 

1.92 

-1.10 

-2.32 

-2.60 

A 1795 

5.63 

2.28 

0.76 

0.08 

A 1068 

5.21 

0.48 

-3.06 

-5.33 

A 1413 

0.65 

-2.42 

-5.41 

-5.66 

A 2390 

3.01 

-0.76 

-3.01 

-4.07 

A 2199 

1.26 

-1.22 

-1.98 

-0.97 

A 2029 

-1.15 

-3.48 

-4.54 

-4.45 

A 2597 

2.86 

-1.17 

-3.82 

-5.44 

Hydra A 

4.07 

1.19 

0.34 

1.05 


presenting a detailed study of three clusters with both 
arcs and X-ray data available, conclude that the mass 
estimates from the arc modelling can be nearly a factor 
~ 2 — 3 larger than those from the X-ray observation in 
the innermost regions. As summarised by Kneib (2000), 
the discrepancy might be due to different reasons: too 
simple X-ray modelling, merging and projection effects. 
Also non-thermal effects can play a role, although mag¬ 
netic fields can be ruled out as error sources (Dolag & 
Schindler 2000). 

Since gravitational lensing estimates provide projected 
masses within projected distances from the centre, the 
X-ray mass must be prcAected too, in order to achieve 
the comparison. In Sect. ^ the X-ray mass was calculated 
within a spherical volume, while for the comparison it has 
to be computed within a cylindrical volume. We investi¬ 
gate projection effects by considering a triaxial /3-model 
and the 4 models used in Sect. H] and |[ As done in Sect. 

we compare the estimates from the triaxial modelling 
with those obtained assuming spherical symmetry, but in¬ 
tegrating Eq. (||) within a cylinder of radius R. We thus 
obtain a projected masses and for the 

spherical model and a triaxial model with a core radius 
along the line of sight c = i x r®^*, respectively. Formally, 
the integral of the total mass density extends from the ob¬ 
server along the line of sight through the cluster infinitely; 
in practice, a cutoff is used. In the following, a cutoff equal 
to 4 X rsoo is used, as a compromise between to much ex¬ 
trapolation and too short integration length. Although the 
integration volume is different from the one used in sec¬ 
tion the ranges of elongations i given in Table ^ apply 
in this case as well. 

Similar to the relative errors defined in Sect. ^ and we 
define: 


p. "/SW - 

^mo\{R) = - 




( 7 ) 


Since mass estimates from the arc modelling are usually 
larger than those from the X-ray observations, a triaxial 
model that implies a negative relative error can allevi¬ 
ate the discrepancy between lensing and X-ray mass es- 



Fig. 5. The relative errors (i?) for the projected mass 
estimates for A2390 plotted for the 4 triaxial models: com¬ 
pressed along the line of sight (z = Zmin) (dotted line), pro¬ 
late {i = e) (dashed-dotted line), oblate (z = 1 ) (dashed 
line) and elongated (z = Zmax) (solid line) shapes. The cut¬ 
off is 4 X rsoQ. Negative relative errors imply underestima¬ 
tions of the total mass if spherical symmetry is assumed: 
this is the case if the cluster is elongated. Overestimations 
are found for compressed clusters. Therefore, an elonga¬ 
tion along the line of sight can contribute to resolve the 
discrepancy between X-ray and lensing mass. 


Table 10. Relative errors for the projected masses at i? = 
j,circ ^ models. The cutoff along the line of sight is 

4 X rsoo. 


Cluster 

r*min 

Proj 

(%) 

(%) 

^proj 

(%) 

ZT'^max 

-^proj 

(%) 

C10016-tl6 

41.54 

11.03 

-3.55 

-41.07 

A 478 

39.43 

8.29 

-8.61 

-40.80 

A 1795 

41.24 

9.71 

-8.81 

-37.81 

A 1068 

40.01 

11.08 

-12.38 

-30.63 

A 1413 

39.10 

14.26 

-16.44 

-20.27 

A 2390 

43.68 

13.70 

-9.01 

-29.71 

A 2199 

38.49 

6.89 

-7.93 

-44.57 

A 2029 

38.75 

9.11 

-11.98 

-34.46 

A 2597 

40.26 

10.47 

-12.24 

-32.12 

Hydra A 

39.68 

6.67 

-7.19 

-46.37 


timates. This is the case if the triaxial model describes 
elongated shapes, as shown in Fig. ||, where we plot the 
relative errors (R) for the cluster A2390 (the trend of 
the relative errors is the same for all the sample clusters). 

In Table we list the relative errors at i? = for 
i = *mirn e, 1 and Zmax- Within this projected distance we 
find the average values: < > 3 = 40%, < E^^^- >s= 

10%, < El^^. >s= -10% and < >s= -36%. In 

Table we list the relative errors at i? = rsoo for the 
4 models. Within R = rsoo we find the average values: 
< >s= 33%, < E;^^. >s= 9%, < El^^. >s= -6% 

and < Ap™ >s= —27%. These results show that X-ray 
estimated projected masses are larger for elongated clus¬ 
ters than for spherical ones and that discrepancies up to 
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Table 11. Relative errors for the projected masses at i? = 
rsoo for the 4 models. The cutoff along the line of sight is 
4 X rsoQ. 


Cluster 

-^proj 

(%) 

(%) 

-c^proj 

(%) 

P^max 

-^proj 

(%) 

C10016-tl6 

40.70 

14.35 

1.65 

-30.89 

A 478 

30.77 

7.31 

-5.70 

-30.72 

A 1795 

29.50 

6.27 

-7.69 

-29.78 

A 1068 

31.44 

9.65 

-8.45 

-22.69 

A 1413 

35.65 

16.04 

-8.78 

-11.90 

A 2390 

31.33 

9.01 

-8.27 

-24.14 

A 2199 

30.42 

6.55 

-4.88 

-33.42 

A 2029 

32.91 

10.20 

-6.35 

-24.18 

A 2597 

32.84 

10.22 

-7.45 

-23.10 

Hydra A 

29.60 

5.04 

-5.51 

-35.57 


~ 30% between the two methods can be resolved if a tri- 
axial /3-model with a maximal elongation /max is used. 
For the cluster C10016-I-16, Small et al. (1995) derive a 
projected mass of 7.3 x IO^^Mq integrated out to a pro¬ 
jected radius of 600 kpc using X-rays, and 8.5 x IO^'^Mq 
from weak leasing. Using the spherical model, our esti¬ 
mate for the projected mass is 6 x IO^'^Mq if no cut-off 
is used and 5.6 x IO^^Mq with a cutoff equal to 4 x rsoo. 
Instead we find that using a triaxial model the discrep¬ 
ancy disappears if the core radius along the line of sight is 
~ 1.76 X (we use a cutoff equal to 4 x rsoo). C10016-I-16 
is among the most X-ray luminous clusters known and X- 
ray selected, thus a shape which is elongated along the line 
of sight is not surprising. From this estimate we conclude 
that elongation is one of the factors contributing to the 
discrepancy between lensing and X-ray mass estimates. 


7. Constraints from the Sunyaev-Zel’dovich effect 


From X-ray observations only, it is impossible to have 
information on the cluster elongation in the direction of 
the line of sight. The uncertainties on the total mass and 
gas mass fraction estimates due to this unknown quan¬ 
tity might be reduced by noticing that the core radius 
along the line of sight can be constrained using a com¬ 
plementary and independent measurement: the SZ effect. 
Assuming the triaxial /3-model as in Sect. (||) one gets: 


y'^{x,y) Hq 
S x{x,y) 




where y is the SZ Compton parameter, Hq is the Hubble 
constant, r®" is the angular major core radius and i is 
defined as in Sect. i) (see Puy et al. (2000) for the 
computation of Sx and y). Thus, the combination of 
X-ray emission and the SZ intensity change, which is 
usually used to constrain the Hubble constant under some 
geometrical assumptions (spherical, oblate or prolate 
shapes), can be used to estimate the elongation along 
the line of sight i, provided that a value for the Hubble 


constant is assumed. Although the SZ effect is strongly 
affected by temperature gradients, we again assume 
isothermality of the ICM because of the lack of detailed 
temperature maps. As an example we estimate i for the 
cluster C10016-I-16 using the measurement by Hughes & 
Birkinshaw (1998). Assuming Hq = 50kms“^ Mpc“^ and 
the value y(0,0) = (2.27 ±0.36) x 10“"*, for the Compton 
parameter measured for the line of sight going through 
the cluster centre we obtain: i = 0.87lg;22- For i = 0.87 
the errors due the assumption of spherical symmetry are 
negligible, since in this case the relative errors defined in 
Sect. (|^) and (||) are less than 3 % for any R considered. 
Taking the errors on the estimated i into account, we find 
the following results. At rsoo, the relative errors for the 
total mass are less than 3% and those for the gas mass 
fraction less than 7%. At r®*’’®, the relative errors for the 
total mass are 8% and —11% (for upper and lower limit 
of i, respectively) and, those for the gas mass fraction 
less than 7%. 

Instead, using Hq = 71 ± 7kms“^ Mpc“^ (Fukugita 
& Hogan 2000), we obtain: i = 1.241 q4o. 

find: ^^^^(rsoo) = 1.7%, EI^^^tqqq) = -1.3%, 

= 10% and = -0.4%. We 

conclude that in this case a triaxial /3-model improves the 
estimation of the total mass within the central region. 
We note that the estimated elongation of i ~ 1.24 agrees 
qualitatively with the conclusion of Sect. but this 
smaller value suggests that elongations along the line of 
sight can only partially resolve the discrepancy between 
lensing and X-ray mass estimates. 


8. Conclusions 

The elliptical X-ray emission shown by many galaxy clus¬ 
ters make us ponder about the possible triaxial distribu¬ 
tion of the ICM. In this paper we model the X-ray surface 
brightness by means of a spherical and a triaxial /3-model 
and, assuming hydrostatic equilibrium and isothermality 
of the intracluster gas, we compare the estimates on the 
total mass and gas mass fraction we obtain from these geo¬ 
metrically different models. Strongly elongated shapes are 
not allowed in this model, suggesting that if these should 
be nonetheless observed they may be indicative of unviri- 
alised clusters. Analysing 10 ROSAT clusters we find: 

— If prolate or oblate shapes are assumed, then the differ¬ 
ences in the estimates of both total mass and gas mass 
fraction are negligible (less than 3%) at every distance 
from the cluster centre. 

— If the cluster is compressed along the line of sight, the 
total mass in the central regions of the cluster is un¬ 
derestimated. The underestimation depends on the de¬ 
gree of compression, the more compressed the shape, 
the larger the underestimation. 

— If the cluster is elongated along the line of sight, the 
total mass in the central regions is overestimated. The 
overestimation depends on the elongation: the more 
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elongated the ICM distribution, the larger the overes¬ 
timation. 

— At large distances from the cluster centre the differ¬ 
ence for the total mass estimates is negligible for com¬ 
pressed and elongated shapes: at rsoo these differences 
are generally less than 4%. 

— The gas mass fractions are overestimated if the cluster 
is compressed along the line of sight. 

— If the cluster is elongated along the line of sight the 
gas mass fractions are slighty underestimated, mainly 
at large radii. 

— We find that projection effects are important when 
comparing X-ray and lensing mass estimates and that 
quite large discrepancies between the two methods can 
at least partially be reconciled if the ICM is elongated 
along the line of sight. 

Since the core radius along the line of sight is not mea¬ 
surable by means of X-ray observations only, we use SZ 
measurements to improve the modelling of the ICM dis¬ 
tribution with a triaxial model. In fact, assuming a value 
for the Hubble constant, we are able to estimate the elon¬ 
gation of the ICM along the line of sight and to constrain 
more precisely total mass and gas mass fraction. In this 
context SZ maps are very important: the modelling of the 
SZ temperature decrement with a /3-model does not only 
allow the comparison of the resulting best ht parameters 
with those from X-ray analysis, it gives better estimates of 
the central Compton parameter than single measurements 
through the cluster centre. This method can clearly also 
be used if the X-ray emission shows circular isophotes. 
Consequently, X-ray and SZ maps together, make the in¬ 
vestigation of selection effects in X-ray clusters possible. 
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